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1. Introduction 

In recent years, considerable progress has been made in understanding plasma 
phenomena in the radiation era (plasma epoch), a period of the expanding universe 
which has particular interest. It is known that after the end of the inflation, 
the universe enters standard Friedman-Robertson- Walker(FRW) stage and the 
wavelength of the large-scale perturbations increases less rapidly than the Hubble 
length which, at that time corresponds to the particle horizon (Zimdahl W 1997). It 
is known that at the end of the inflationary era density perturbations and gravitational 
waves(GW) are produced (Brevik I and Sandvik H B 2000, Lyth D H 1996). The 
density perturbations are thought to be responsible for the formation of large-scale 
structure along with a possible GW contribution. Moreover, the density perturbations 
due to the inflation are responsible for the observed anisotropies of the Cosmic 
Microwave Background Radiation(CMBR) (Kandaswamy Subramanian and Barrow J 
D 1998, Koh S and Lee C H 2000). On the other hand, it is argued (see Dimopoulos 
K et al 2001) that at the end of the inflationary epoch large-scale magnetic fields 
and inhomogeneities could be generated. The existence of large-scale magnetic fields 
in the early universe may have important consequences on the large scale structure 
formation process as shown in (Maartens R et all C 2001, Tsagas G C 2001, Grasso 
D and Rubinstein H R 2001, Tsagas G C and Maartens R 2000). The history of 
the magnetic fields from magnetogenesis in the early universe to the present time 
with the evolution and damping of the magnetic fields up to recombination era, is 
presented in (Jedamzik K et al 2000, Olinto V 1998). All this knowledge related 
to the density perturbations and magnetic inhomogeneities in the early stages of 
the universe, help us to construct possible scenarios to discuss plasma phenomena 
(see Holcomb K A and Tajima T 1989, Holcomb K A 1990, Dettmann et al 1993, 
Sahni V 1989, Ford L H and Parker L 1977, Caprini C and Durrer R 2001, Brodin 
G, Marklund M 1999, Marklund M et all 2000, Brodin G et al 2001, Brodin G et 
al 2001). In particular, from t = 10~ 2 s to the time of recombination the ordinary 
matter content will be in the form of plasma. From approximately t = 10~ 2 s to 
t — Is, the plasma is dominated by electrons and positrons while from t — Is to the 
time of recombination at roughly t = 10 13 s, the plasma consists mostly of electrons 
and protons with a admixture of ions of light elements. These plasmas are believed 
to be in thermal equilibrium, or at last strongly coupled, with photons. It is the 
photons however, which dominates the energy density of the universe until the time 
?s 10 2 s Thus most of the period prior to the recombination is traditionally called the 
"radiation era" of the universe. In this paper we will focus on the radiation epoch of 
the universe and try to discuss, whether the existence of the GW in a FRW space time 
with zero spatial curvature and small magnetic field may produce instabilities. These 
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kind of problems are of particular interest since astrophysicists have began to discuss 
finite-amplitude wave propagation in hydromagnetic systems in the presence of GW 
(Papadopoulos D et al 2001, Papadopoulos D Vlahos L and Esposito F P 2001). 
Long ago Papadopolos and Esposito (Papadopoulos D Esposito F P 1982), derived 
the exact equations governing finite amplitude wave propagation in hydromagnetic 
media within the general theory of relativity, in the so called Cowling approximation 
(Cowling T G 1941). These equations are of particular interest because when they 
are linearized, they are very useful to discuss stability criteria in the presence of 
magnetic fields. Recently, (Papadopoulos et al. 2001) have derived a dispersion 
relation in an anisotropic cosmological model (Thorne K S 1967). They verified 
that in the case that the parameter 7, used to parametrize the equation of state, 
approaches 1/3, the anisotropic cosmological model becomes isotropic, actually a 
FRW with a scale factor S(t) = t 1/2 , and the expansion 6 of the universe prevents 
any instability. But in ((Dimopoulos K et al 2001) the authors argue that in FRW 
cosmological model, if one starts with small fluctuations, may obtain growth in those 
instabilities. The classical theory of GW perturbations on FRW universe has been 
worked out by Lifshitz and discussed by several authors (Ya.B Zel'dovich and Novikov 
I D 1983 ). In this paper we consider a FRW cosmological model,with signature 
(— , +, +, +) witch is spatial flat and a weak magnetic field, perturbed by a GW with 
two polarization modes. As state by Marklund M et al 2000, the presence of a weak 
magnetic field propagating on the background FRW model couple to the gravitational 
wave producing a pulse of gravitationally induced electromagnetic wave. In principle, 
this breaks the homogeneity and isotropy of the FRW spacetime. However, it can be 
shown that a weak magnetic field is consistent with such a background provided that 
the spatial variations it generates are negligible. Therefore, we assume that both, the 
weak magnetic field and the GW, do not break the homogeneity and isotropy of the 
FRW space time. Applying the formalism presented in (Ya.B Zel'dovich and Novikov 
I D 1983, Brandenburg A et al 2000, Weinberg S 1972) in our scenario, we verify that 
density fluctuations may be obtained. 



2. The Equations 

The exact equations governing finite-amplitude wave propagation in hydromagnetic 
media in the frame of general relativity have been discussed in (Papadopoulos D 
Esposito F P 1982, Papadopoulos D. et al.2001, Papadopoulos D Vlahos L and 
Esposito F P 2001). For completeness we recast the relevant equations. We start 
with the Einstein field equations 

Rab — -^9abR = —KT ab , (1) 
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Tf = (2) 



For simplicity we obtain c = 1 and k = = 1. Taking the covariant divergence of 
the Bianchi identities we obtain 

{R ab -\Rg ab U = u (3) 

where R ab = R c abc is the Ricci tensor and R a bcd is the curvature tensor defined by 

Ua-bc ~ U a . cd = U d R d ahc (4) 

For a unit time-like vector we choose u a u a = — 1 and our hydromagnetic system will 
be specified by the following choice for the energy-momentum tensor 

T ab = (e + ^-)u a u b + {p + ^-)h ah - H a H b (5) 

with 

h ab = g ab + u a u\ e = p + pU (6) 

where u a is the fluid velocity, p the mass density, pIT the internal energy density, p 
the pressure of the fluid and H a is the prevailing magnetic field as measured by an 
observer co-moving with u a . Furthermore we introduce the expansion 9 = u a a , the 
shear a a b = h a hfu^ c; d) — \0h a b and the twist u a b = ha^b u [c;d], where the round bracket 
denotes symmetrization while the square bracket antisymmetrization. From Eq.(l) 
and Eq.(5) we find 

R ab = - {he + 3p + H 2 )u a u b + he-p + H 2 )h ab - H a H b } (7) 
Substituting Eq.(7) into Eq.(3) we have 

x + g ab (p + ^P) + 2*9 + x, a ii a + x(9 + 9 2 + u%) - (H a H b ) ;ab = (B) 

V 1 J -ab 

where ii a = u a c u c and x = e + p + H 2 (Notice, that Eq. (8) can be obtained directly 
from Eq.(2) and Eq.(5)). Raychaudhuri's equation in the form 

u% = e + j + 2(a 2 -uj 2 ) + 1 -(e + 3p + H 2 ) (9) 

where 2a 2 = cr ab cr ab , 2uo 2 = uo ab uo ab and dot means covariant derivative along u a allows 
us to write the Eq. (8) as 

x + g ab (p + — \ + 2x9 + 2x9 + x. a ii a 

V 2 / -ab 

+2x(— + a 2 - u 2 ) + -x(e + 3p + H 2 ) - (H a H b ). ab = (10) 
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In order to make further progress with Eq.(lO) we must use the equations of motion 
for the fluid and Maxwell's equations. These equations are 

tt2 

= Tf = ±u a + xii a + x6u a + (p + — ). b g ab - (H a H b ), b (11) 

and 

H a = (a a b + co a b - \5 a b 6)H b + -^-p ;b H b u a (12) 
o a -\- p 

The last equation (12) may be written in the following form 

where fi is the permeability. Equations (11) and (12) imply that 

e + (e + p)6 + (e + p)6 = (14) 
Recalling that x = e + p + H 2 we find that Eq.(lO) takes the form 

(e - ^-)-,a b u a u b = h ab (p + ^-) ;ab + 2(H*e) - (H a H% b 

29 2 x 
+2x(— + a 2 - uo 2 - ii a u a ) + -(e + 3p + H 2 ) 

+2u a {H a H b ), b + {H 2 ), a u a (15) 

The Eqs. (11), (12) and (15) may be applied to the investigation of perturbation 
effects and hence to the study of the linearizied stability criteria. 



3. The Perturbed Equations 

Following the well established method by (Hawking 1966), used also in (Papadopoulos 
D Esposito F P 1982, Papadopoulos D et al. 2001, Papadopoulos D Vlahos L and 
Esposito F P 2001), we consider the perturbations as follows: 



e = e + Se,p = p + Sp (16) 

v? + = H% + 5H» (17) 

where the index means unperturbed quantities. 

Taking into account the Eqs. (16), (17), we perturbed the equations Eqs. (15), (11) 
and (12) respectively, keeping only linear terms in the perturbed equations and 
assuming that 5g ab ^ 0. Thus we have 

5[(e-^), ab }u a u b -5T k ab (e-^-), k u a u b 
+ (e - \)Au a tu b + u b 5u a ) + (Se - H c 5H c ) ]ab u a u b 
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= 5g a % h + (u a 5u b + u b 5u a )p* ab + h ab {5p* ah - 5T k ab p* k ) 

+ 25[(H 2 ), C Q + H 2 Q, c ]u c - H b [28{H« ab ) + H k 5{T a aKb ) + 28V a ak H k b \ 

- 2H a ab 5H b - 25T a al H l H b b - 5(H a b )H b a - Htf(H b a 

- 5Y a bl H l H b a - H a H c 6T b ^ b - 5T b ac HlH a + 2x[a5a - u c u d {5u a c u a4 + u%6u a>d 
+ 5Y a d u l u a4 - u%u m 5TZ) - (u d 5u c + M C <5«» M ] 

+ -{e + 3p + H 2 ){5e + 5p + 5H 2 ) + -x{5e + 35p + 5H 2 ) 

+ [25« c u 0i c + 2u c (5u a , c - 5T2u m )}(H a 5H b + H b 5H a ) + (H 2 ). a (5u a ) ;c u c 
+ (H 2 )XM C + (H 2 ) ta u c (6u% + 5r^ Mm ) (18) 
The perturbed Maxwell's equations are 

8H b + 5T b ol H l = - 5a b H a + ^H b 56 - ^Sp, c H c (19) 
and the perturbed Eqs.(ll) give 



- (e + p + H 2 )5u c fi = h c a [H b 5H% + 5T a bl H b H l 

+ H a 5H b b + 5T b bm H m H a ] - 5h ch p* b - h cb 5p* b (20) 

We will apply these equations to a spacial flat FRW space-time given by the 
equation (see Ya.B Zel'dovich and Novikov I D 1983) 

rfs 2 = _ c 2 dt 2 + 5 2 [(1 + ^ )(fe 2 + (1 _ h )d 2 y + ^] + 2/^^y ( 2 1) 

where S = S(t) is the scale factor of the FRW universe, parametrazing the universe 
expansion with hi = h+(t,z),fi2 = h x (t,z) are the two components of the GW 
corresponding to the two possible polarizations. For the metric (21) the non-zero T c ab 
are 

o 1 o 1 



1 

2 

„1 5* o 1 / hi n hiS n 

= + -f-i# - 2^^ 



5 2 l S 2 S 3 h 



1 02 - 52 z 5<3 ^ 

r i _ 1 fc r i _ 1 fc p2 _ n h 2 S fi 

1 13 - 2^52 ftl ^ ' 1 23 - 7^2 n 2,z0, J- 01 - + o" ^S 2 " ~~ g3 

F 2 _ ^0 _ 1 ^1,0 _ 9 ^1^0\ F 2 _ p 2 

~~ S 2 S 2 S 3 ~~ 2 S* 2 3 ~~ 2S 2 

F 3 _ ^0 p 3 J_, p 3 j_, p 3 _J_, / 99 x 

03 S n ~ 2S 2 2 ~ 2S 2 z 22 ~ 2S 2 z 
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In an expanding universe described by the metric (21), we obtain S(t) = t 1 ! 2 and 
consider an equation of state p = |e(radiation era). We assume that u a = (— 1, 0, 0, 0), 
H a = j5(0,H%,Hg,H$) (Holcomb 1990), with H% = const,H 2 = (Hi) 2 + (H 2 ) 2 + 
(Hi) 2 , 5H a = jz(0,5H\5H 2 ,5H 3 ), 5u a = (5u° , 5u l , 5u 2 , 5u 3 ) and k a = (0,0,0, A;) 
where k is the magnitude of the comoving wavelength (momentum). Thus, the 
perturbed Maxwell's equations (19) take the form 

6H* = 9SH 1 + H l 50 - heH 1 - H] Q ))8u° 

5H 3 = 65H 3 + H 3 56 - ^9H 3 5u° - -if MuJ, - H 3 5u% (23) 
and the equations of motion (20), reduce to the equations 
6u° = 



x6v} = \x65u 1 - H 3 5H\ - ffWg - ^jf ^ h 



1,3? 



x5u 2 = \x65u 2 - H 3 5H 2 3 - ^\ 2 , 3 , 



2 -x65u 3 - 2H 3 5H% + i 
In this case the Eq.(18) becomes: 



x5u% = -x65u 3 - 2H 3 8H% + — 2 H 0li 5H% + 6p, 3 + (H v ) 2 h^ (24) 



3 



-OSp,o + 3—^ [^g + 2(e + 2p + — )} ^ ) 

-4056(e + P~§- 2 ) + ff^,o - ^(5e + 5p) 

~(Se + Sp)(e + 3p) ~\{5e + 35p)(e + p)-^(Se + 35p) 

= (#TM-§ + e + 2p+ff-^] (25) 

It is known that in an expanding universe the perturbations can not be solutions for 
the GW, because the amplitude of such a wave must decay in time and its frequency 
must redshift. Thus, writing the wave equation for a free wave in a cosmological 
background and requiring a propagating solution, the spatial part has the usual 
form exp (ikz), while the time depended part has the form exp (iut/S). In the 
plasma period, we assume that the perturbations, Eqs.(16),(17), have the form (Ya.B 
Zel'dovich and Novikov I D 1983): 5H a = exp (ikz — iujt), 5e = ^4 exp (ikz — iuit), 
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5u a = 5uq exp (ikz — iut) and the two components of the GW (Brandenburg A et al 
2000, Weinberg S 1972) hi, h 2 = ^exp (ik g z — iuj g t) where uj = ^, tu g = ^ and c^, 
ujgi are the frequencies at the initial time t; L < t. Under these considerations, Eq.(25) 
yields a general dispersion relation for the coupling of GW with MHD waves which 
is: 

AR^e + ±-H 0a 5H a R 2 = h^H 1 ) 2 ^ (26) 



where 



A^U-^) + f(l + 4 " ; 



2 V g 4 > 1 gV ' S 4 I 

cu 2 . u 2 A . ico „ 2 ,u 2 A . k 2 c 2 

- l -{e + 3p)(l + c^) - i(e + p)(l + 3c 2 ) - + 3c 2 ) 

uj 2 . u 2 A . 8iuu 2 A k 2 Hi. n2 / 2&u A 47 ^ 

«" = T (1 -#— wr + ^-^ + ^ + ^ + ^T^-T' 

i?3 = -^ + e + 2p+-^-^ + 8^ 2 
«A f w 9 2 9^ 16 2 

and = -rrV, 1/4 = — is the Alfven speed and c s is the sound speed. It is important 
to emphasize that, the second component of GW does not appear in Eq.(26). Besides, 
the right hand side of Eq.(26) is proportional of H 1 . This suggests that the dispersion 
relation (26) may be used for waves that correspond to the magnetosonic mode 
(H 1 ^0,tf 3 = 0). 



4. The dispersion relation 

Furthermore, obtaining H l ^ and H 3 = 0, the Maxwell's equations and the 
equations of motion give the following results respectively: 5H 3 = 5H 2 = 0, 

WSH> = - f )d + cj)fc + - f ) + f ] (28) 

and 5m 1 = 0,<5w 2 = 0. We notice that in all equations above equations we have taken 
c = 1, k = %%Gjc SL = 1 and we have put H = B — > \^AttH (from Gaussian units), 
where n is the coupling constant in Einstein's field equations. Substituting Eq.(28) 
into Eq.(26) and using ordinary units, we obtain 

( H 1 ) 2 

(A r + iAi)5e = hy^—!-{B r + iB t ) (29) 
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where 

A- = -(^) 2 s 2 + S , A = -(^) 3 S 3 + (^)Si, 

fi r = -(|) 2 n 2 + ^n 1 + n , 



with 



and 



B t = -{-yM 2 + (-)nM l + M (30) 



—3 



c 2 u\ 2u 2 A u\ c 2 u\ 



--d-H)[^-^ + 3 P )( 1 + |) 

_^ (£ + p)(1 + 3 £i)-^ (1 + 3|), 

+351(1 + -f )[^ - — ( £ + 2p+ ^(1 + -§)(! - 25^)] 
So^«l + ^)[^-^(e + 3 rt(1 + |) 
^ (e + p)(1 + 3 |)-^( 1 + 34)l 



c 4 v c 27 c 4 S 4K c 2 ' 

u\ c 2 k 2 c 2 IQitG. H 2 q2 ,28u 2 a 47 



(l + i)[^ 2 ---^(^ + 2p + ^i) + ^(^-^)]} (31) 



c 2 S 4V c 2yL S 2 c 4 v ^ 4ttS 4 ' v 9c 2 S 4 3 

112 " g^ 1 -^)^-^) 
ni = 50[(l-^|)(l-|-^ 8 ^ 



c 2 5 4yv 5c 2 5 4 ' 5c 2 5 4J 
r/ „ Au 2 A , r k 2 c 2 u 2 A uo 2 87TG, n H 2 



gLV s2 4c 2£ 4 c 4 V 4?r>s4 2 

lfaG ( e + 2P+ ^) + ^(H|-?)]} (32) 



c 4 v v AttS^ v 9c 2 S 4 3 



M 2 = ^(l U 

2 0„ 2 04V 



2c 2 S 4> c 2 S 4 ' 
„ 2ti -2 ,, fe 2 »U StiG, „ ff 2 

= (i - ^[-ir - + — ( f + 2 p + ^) 
- f|)] 
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From now on, we will use the notation n = ^ and n g = ^j-. 

It is evident that this section 4, from the beginning, relies heavily on equations 
and calculations, therefore we shall try to present our case by splitting this section 
in two subsections, one dealing with the gravitational-wave-free case and another one 
addressing the impact of the gravitational wave. Notice that the coefficient of the 5e 
on the left hand side of Eq.(29) describes the dispersion relation of the fluid. While, 
the coefficient of hi on the right hand side of Eq.(29) describes the dispersion relation 
of the gravitational wave. 

4a. The gravitational-wave-free case 

In the absence of the gravitational wave {hi = 0), the dispersion relation 
(Eq.(29)) is: 

A r + iAi = (-n 2 S 2 + S ) + i(-n 3 E 3 + nE ± ) = (34) 

Eq.(34) has three roots, two complex conjugates and one real; we examine one of 
them, say the rii = n r + irii, with rii « n r , since its complex conjugate will give 
the same information regarding the fluid and the real one is not important from the 
physical point of view. We substitute ni in Eq.(34) and the resulting equation after 
neglecting non-linear terms of is: 

- {n 2 + 2m i 72 T )S 2 + S - i{nf + 3mjn^)S 3 + i(n r + in^Ei ~ (35) 

From the imaginary part of Eq.(35), we obtain n 2 ~ — 2n^S| + S^-, and rii ~ ^^, 2 ~^^ . 
After some elementary calculations, we verify that the condition n« << n r is satisfied. 
Thus, from Eq.(29) in the absence of the GW, we obtain Eq.(34) and find the root 
rii which has a real part n r , given by the equation 

2 _ 81 i- 3 -g5-J C s , * K U A 

r ~ 4. 2^4 , 40091] S 2 3S 6 { J 

I S 2 T 3 S« ^ 81 J 

and an imaginary part, given by the equation 

n t ~ s 2 + I gg \ (37) 

,k 2 c 2 , 2 k 2 u\ 40002 1 V I 

V S 2 ~T 3 S 6 ^ 81 J 

Note, that the two frequencies, n r and rii, depend on the expansion 9. If, somehow, the 
universe would stop to expand, then the real frequency n r is just the so called magneto- 
sound frequency (Jedamzik K et all 2000, Olinto V 1998), while the frequency rii 
becomes zero. In an expanding universe, an imaginary part of the frequency rii 
appear, and the real frequency n r , is shifted by a term proportional to 9 2 . The 
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existence of the frequency n r means that before the GW starts to interact with the 
plasma, the plasma were oscillating with this n r normal frequency, while the frequency 
rii corresponds to just a negligible noise. 

4b. The dispersion relation in the presence of the gravitational wave. 

Now we consider the case where the GW interacts with plasma. We assume 
that the GW is weak and does not affect the background metric and the frequencies 
n r and rij. In this scenario, we write: 5e ~ -^Seoexpi(nt — kz), and hi ~ 
^hio expi(n g t — k g z) (Brandenburg A et all 2000, Weinberg S 1972), where 5eo = 
constant and h w = constant. Thus, Eq.(29) is written as 

5e ~ fc^(^M).l|[^i] (38) 

The term T = J±gi takes its maximum value close to the frequency n — n\. Let 
n g and k g , the frequency and the wave number of the GW respectively. Close to the 
frequency n±, we assume that the frequency of the driving GW n g coincides with the 
frequency n r = Re(ni). The produced magnetosotic wave has n r = n g = k g c and the 
two wave numbers differ by an amount Ak, e.g. k = k g + Ak. Taking into account 
all the above considerations, Eq.(38) reads 

~ h r -^Ak njt Hi r B g B k 

6e °~ h " e 6 ^AkA k + A-J (39) 
In the linear theory of perturbations we have 8e < e . Thus, Eq.(39) gives 

^ioe- lzA V^[-^-] < e = const., as Ak - 0. (40) 
47TD' 3 AkAk 

B 9 = ^Un) - 8^LA) - ni (-^_ + JUL) 

■ r 130 , k 2 c 2 3k 2 u 2 A 50\, 

+ i[— n t n g + n g {-^ + - — )] (41) 

5* = |[^-^]+<n fl (^) (42) 

^ = ~ n ^ + ^) + ^(^T + ^) ( 4 3) 

Because of the form of (Eq.(37)) the quantity e*™ 1 approaches to a constant value, 
as t approaches to very large values. Thus, to the approximations we made, the 
disturbance 5e does not exhibit growth. However, the term e trii enhanced the 
amplitude of the oscillation and this may cause further enhancement on a non-linear 
approach which we will not discuss in this work. (The result is in agreement with the 
result obtained in (Holcomb K A and Tajima T 1989, Holcomb K A 1990). Alfven or 



where 



and 
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other MHD oscillations are not enhanced, but magnetosonic frequencies through the 
interaction with the GW become dominant. Therefore, it is interesting to know the 
maximum length which the GW will travel in the plasma and the maximum time the 
GW will interact with the plasma. From Eq.(40) we have 



»S <Re(B s A- t ) 



Eq.(44) can give us approximately the Ak and subsequently, we can calculate the 
maximum distance L over which the two waves may interact coherently. That 
maximum L can be calculated in terms of the wavelength \ g through the relation 
L = A„(l - f ) e.g 

£<A,{l-wffi ^fl - ]} (45) 



where 



3S 2 v 9 S 6 ' L S 2 v S 2 S 6 
^.tf,-^^)] (46) 



and 



2k g c 2 Ak g u\ IWn, k 2 c 2 3k 2 u 2 A 56 2 
h - (— + -35r)(— - -gT + -g* r ) (47) 

/3 = ^-2n l( ^S + ^) (48) 
( 2k g c 2 s 4:k g u 2 A 2 

/4 - + (49) 

The right hand side of Eq.(45) is smaller than the particle (Peacock J A 2000) horizon 
of comoving radius rn — c j Hence, in a small fraction of the universe radius, there 
exists an interaction between the GW and the plasma producing small fluctuations 
with amplitude P = h 10 exp i[nit — A/^]-g§-[ ^ r ^' ]. During the interaction, energy 
of the order P 2 L 3 is transferred from the GW to the fluid and magnetosonic waves 
are produced. These magnetosonic waves do not get damped since in our model we 
do not include dissipating processes. In a future paper we intend to study a more 
realistic model and investigate further astrophysical implications. 
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3. Discussion 

In the present paper, we have applied the equations governing finite amplitude 
wave propagation in hydromagnetic media discussed by Papadopoulos and Esposito 
(Papadopoulos D Esposito F P 1982), in the so called Cowling approximation in 
a perturbed spatial flat FRW cosmological model. We have assumed that, despite 
the presence of the magnetic field and the perturbation h±, the homogeneity and 
isotropy of the considered cosmological model do not change. Upon the consideration 
of those equations, we have derived a dispersion relation (Eq.(29)) in the radiation 
era (plasma epoch) . We found that, if the gravitational waves and the magnetic field 
(produced after inflation period) are parallel, then the gravitational waves may excite 
fluctuations which appear as fast magnetosonic waves. From Eq.(29), it is evident 
that both, the magnetic field and the gravitational wave should be different from zero 
otherwise we do not have genaration of density fluctuations. The amplitude of the 
magnetosonic wave is proportional to the gravitational wave and does not exhibit 
a damping since in our model we do not take into account dissipative processes. 
The process we suggest seems to be reversible in the sense that energy is transfered 
from the gravitational wave to the plasma exciting magnetosonic waves and vise 
versa. The transfer of this energy from one side to the other occurs without loses 
and is of the order P 2 L 3 . Also, we have found that in an expanding universe, in 
the absence of the gravitational wave, a new imaginary part of the frequency (rij) 
appears and the real part is shifted by a term proportional to 9 2 . In the presence 
of the gravitational wave, the interaction of the plasma with the gravitational wave 
produces magnetosonic waves oscillating at a frequency n r = n g = k g c. Finally, 
the application of the above equations to the study of small disturbances is rather 
straightforward although tedious at this time, but we hope to be able to discuss, in a 
future publication, these equations in more realistic hydromagnetic media. 
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